Heat exchanges in a quenched ferromagnet 
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The off-equilibrium probability distribution of the heat exchanged by a ferromag- 
net in a time interval after a quench below the critical point is calculated analyti- 
cally in the large- iV limit. The distribution is characterized by a singular threshold 
Qc < 0, below which a macroscopic fraction of heat is released by the k = 
Fourier component of the order parameter. The mathematical structure producing 
this phenomenon is the same responsible of the order parameter condensation in the 
equilibrium low temperature phase. The heat exchanged by the individual Fourier 
modes follows a non trivial pattern, with the unstable modes at small wave vectors 
warming up the modes around a characteristic finite wave vector kw Two internal 
temperatures, associated to the k = and k = ku modes, rule the heat currents 
through a fluctuation relation similar to the one for stationary systems in contact 
with two thermal reservoirs. 

PACS numbers: 05.40.-a, 05.70.Ln, 05.70.-a 

Finding the principles underlying the probability measures of non-equilibrium fluctuations is 
one of the most challenging and far reaching open questions in modern statistical physics. Large 



deviation theory, as recognized in the last years, gives a general theoretica 



scribing probability distribution functions (PDF) in non-equilibrium states However, in spite 
of recent important developments g expneit calculations, especially for interacting systems, are 
limited to few cases A major advance for non-equilibrium stationary states has been the 
recognition of a general symmetry of the PDF of certain observables, described by the so-called 



framework for de- 



fluctuation theorems [3 
temperatures (3% < fa [5 



For example, for a system in contact with reservoirs at two inverse 
6j, the probability distribution P(Q) that the heat Q flows from the first 



2 



to the second heat bath in a large time interval is related to P(—Q) by 

= -AW- (i) 

Fluctuation behavior in non-stationary states is by far less understood. In the thoroughly 
investigated field of aging systems, such as quenched ferromagnets or binary mixtures, disordered 
materia! 

a 



models 



s^and glasses, heat PDF have been considered only numerically in some specific disordered 
and, recently, in an experiment for a brownian particle jg] trapped in an aging bath 9|. 



Understanding the properties of heat fluxes in these systems is of great importance also for what 
concerns the notion of an effective temperature [k]], which is expected to regulate such fluxes 
similarly to what the ordinary temperature does in equilibrium. 

In this Letter we address the latter category of problems, by studying the probability distri- 
bution of heat exchanges in a ferromagnetic model quenched from a disordered state to a final 
temperature below the critical point. We do this by an exact calculation carried out on the time 
dependent Ginzburg-Landau model with an ^-component vector order parameter in the large- iV 
limit H]. Specifically, we find the analytical form of the probability distribution P(Q, t, t w ), where 
Q is the heat exchanged during the time interval [t w ,t] following the quench. Most interesting is 
the existence of a singular threshold Qc, such that for Q < Qc the macroscopic amount of heat 
Q — Qc is entirely released by the zero wave vector mode. This comes about through the same 
mechanism responsible of the transition to the low temperature phase in the equilibrium version 



of the model 
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131 ] . Furthermore, we find that P(Q, t, t w ) asymptotically obeys a fluctuation re- 



lation akin to Eq. ([T]), even though the system is not in a stationary state. This can be interpreted 
as due to the heat exchanged between the condensing k — mode, lowering the system energy as 
an effect of the ordering process, and the modes at some finite characteristic wave vector ku- The 
two inverse temperatures playing the role of /3i,/?2 in Eq. (Tj[|) arise as the typical energy scales 
associated to these two kinds of non-equilibrium modes, and reduce to the bath temperature when 
the system is in equilibrium. 

We consider a system of volume V, described by the Ginzburg-Landau Hamiltonian 



H[<p]= / d d x 



(2) 



'v 

where r < 0, g > 0, and ip = (<pi, ..,<pn), is the iV-component order parameter field. Dynamics is 
governed by the Langevin equation 



dip a /dt = -5H[cp]/5(p a + rj a , (3) 

where <p a ,r] a stand for ip a (x,t),r] a (x,t) and the latter one represents the Gaussian white noise 
generated by the thermal bath with averages < rj a (x,t) >= and < r) a (x,t)r]p(x' ,t') >= 



3 



2T5 a p5(t — t')8{x — x'). The leading order of all quantities of interest in the 1/iV-expansion 
can be obtained by replacing the above Hamiltonian with the time-dependent effective one 



H[tp] 



d d x [(V^) 2 + (r + gS(t))cp 2 } - (NVg/A)S 2 (t) 



(4) 



where S(t) = (ip 2 (x,t))/N must be computed self-consistently and angular brackets stand for the 
average over both initial condition and thermal noise. Due to space homogeneity, the quantity S 
only depends on time. The remarkable feature of the large- iV limit is that the dynamics generated 
by "H, although retaining all the relevant features of the phase-ordering process, becomes exactly 



soluble 
of motion 



lit 



By Fourier transformation one obtains a decoupled set of formally linear equations 




Of 



<p a (k, t) = -u(k, t)<p a (k, t) + Tj a (k, t) 



(5) 



where the stiffness of each mode is given by u(k,t) = k 2 + r + gS(t) and < rj a (k,t) >= 0, 
< r) a (k,t)r)p(k' ,tf) >= 2T6 a p5(k + k')5(t — if). Integrating Eq. (J3J) and taking averages, the 



various observables can be obtained 



151 ] . In particular, the two-times structure factor C(k, t, t v 



(l/V)((p a (k,t)ip a (—k,t w )) with t w <t will play a relevant role in the following. 

The probability to release the heat Q per component in the time interval [t w ,t] is defined by 

KM , H[<p w ], 



p{q- t w , t)= / d[<p t }d[<pJn<Pt, vMQ 



(6) 



N N 11 

[\ k a Vk, a (fa(k, t), <p a (k, t w )) is the joint probability of the two configurations 



where V(<p t , <p w ) 
(ip t , ip w ) at the times t and t w . For a Gaussian process 

C tip WW 



V kja {(Pa,t, <Pa,w) =M" 1 exp 



<pj,t 
a, 



in 2 
ra,w 

^ mil) 



CttC WW 



(7) 



2V(C tt C w 

where the short notation C w t = C(k, t, t w ) (and similarly for Cu, C ww ) has been introduced and M 
is the normalization. Next, using the representation S(x) = e~ zx dz of the Dirac 5-function 

and carrying out the integration in Eq. ()6]), we find 



Zg—lOO 



dz 
27ri 



D Vh{Q,z;t,t w ) 



(8) 



where Q = Q is the heat density and the real quantity Zq is chosen in such a way that the integral 
is well defined 12|. For large V discrete sums over wave vectors can be replaced by integrals, 
yielding 

idi 



h(Q; z\ t, t w ) = -z[Q + gU (t, t w )} - (1/2) J ln[l - zq(k, t, t v 



z 2 b(k,t,t w )} (9) 



4 



where 

U{t,t w ) = [S 2 {t)-S 2 (t w )]/4, (10) 

b(k,t,t w ) = u(k,t)u(k,t w )[C(k,t,t)C(k,t w ,t w ) - C 2 (k,t,t w )}, (11) 

q(k, t, t w ) = u(k, t)C(k, t, t) - u(k, t w )C(k, t w , t w ), (12) 

and the symbol j A denotes an integral with an ultraviolet cut-off A related to the lattice spacing. 

Eqs. ([SEP are completely general. Different dynamical protocols are encoded into the cor- 
relation C. We start our analysis from the simpler case in which the system is in equilibrium 
at a generic temperature T. In this case u(k,t) = u eq (k) and C(k,t,t) = C eq (k) do not de- 
pend on time due to stationarity, so from Eq. ( Tl2l) q = and similarly U = 0. Moreover 



C(k,t,t w ) = C eg (k,t- t w ) = (T/u) eq )exp[-U) eg (t - t w )] |lj|, hence b = T 2 [l - exp[-2u eq (t - t w )] 



In the large V limit, the integral in (JSJ) can be computed by the steepest descent method. The 
saddle point equation , = reads 



d d k z*(t, t w )b(k, t, t w ) , 
(2ir) d l- z *( t ,t w ) 2 b(k,t,t w y 



A 



In order for h in Eq. fl9]) to be defined, it must be z ^ b{k,t ) t w )~ 1 ^ 2 Vk,t,t w . This, for t — t w 3> 
l/weq(A), requires —/3 < zq < /3, with = 1/T. From the above expression of b one sees that the 
integral approaches infinity as z* — > ±/3 so that Eq. ffT3|) admits a real solution z* for any value of 
Q. The large deviation function defined by 

P(Q; t, t w ) = exp{VC{Q; t, t w )} (14) 

is given by C(Q; t, t w ) = h(Q; z*;t — t w ) and is plotted in the inset of Fig. [TJ It is symmetric and, 
for Q not too small behaves linearly, £ ~ — /3\Q\, since z* rapidly converges to (3 as Q increases. 
Notice that the equilibrium temperature can be read out from the singular points of h, which in 
turn regulate the exponential decay of the tails of P. 

Next, we consider the quench from infinite temperature to T < T c , starting with T = 0. 
Since q(k, t, t w ) will now play a central role, let us comment on its physical meaning. Using 
the normal modes decomposition, the average energy at the time t can be written as ((H) = 
f Efe^(M)C(M,*) - (NVg/A)S 2 (t). This shows that q in Eq. flT2J can be interpreted as the 
average heat (per component) exchanged by the individual modes, since the contributions due to 
VgS 2 /4: become negligible at large times as we will show below. 

In a zero temperature quench, one has C(k,t,t w ) = y/C{k^t, t)C(k, t w , t w ), which implies 



6 = 0. In what follows we will consider the large t w limit. In this regimes one finds 
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14 



15| 
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S(t) = —r/g — d/{Agi) and the dynamical scaling property L(t w )~ d C(k,t,t w ) = C(x,y), where 
x = t/t w , y = ktw 2 , L(t w ) = (2t w ) 1 / 2 is the characteristic lengthscale at the age t w of the system, 
and C(x,y) = (47r) d / 2 (— r/g)x d l A exp[— y 2 (x + 1)]. Hence, using these results, also q can be written 
in the scaling form 

= {xy i _ * )x l-i e -*rt - (y 2 - -\e- 2 ^ (15) 

Qtyp 4 4 

where qt yp (t w ) = — [r(87r) ^tt! 2 1 / g] is the typical age-dependent scale of heat fluxes. Notice that, 
for fixed x and y, q grows like tw -1 . Therefore, using the expression of S given above, the extra 
term (NVg/A)[S 2 (t) — S 2 (t w )} oc (t" 1 — t" 1 ) oc t^ 1 (l — x~ l ) is negligible with respect to q, as 
anticipated. The quantity q is plotted in Fig. [2] against y for two different values of x. For any x 
there is a negative minimum —0Q 1 = q(x, 0) at the origin. Since q is the average heat exchanged 
by the single modes, this means that the components around k — cool as the time goes on and 
that the cooling increases with the time difference, as intuitively expected. However, the shape 
of the curves shows that the rate of cooling decreases as y increases, with the unexpected and 
quite interesting feature of the development of a positive peak f3 M l = q(x,yM), which is more 
pronounced for the larger time differences. This implies that the modes under the positive peak 
warm up as time goes on. Since the thermal bath is at zero temperature, this extra heat can 
only originate in the heat redistribution due to the coupling among the modes. In fact, it should 
be kept in mind that the linearization of the equations of motion is only formal, the nonlinearity 
having been preserved through the mean field term S(t). For yet larger values of y the curves 
become flat about zero, indicating that the large k modes are equilibrated. 

In order to see what are the implications of the above features on the properties of P(Q; t, t w ), 
let us compute the integral in (jS]). Recalling that 6 = 0, for h in Eq. Q to be defined, it must be 
z ^ q(k, t, t w )~ x V7c, t, t w . With the form ([15]) of q (see Fig. [2]) this translates into — f3 < z < (3 M . 
The analyticity domain of h is shown in the inset of Fig. [2j Notice that, in this far from equilibrium 
situation, /3q and (3m play a role analogous to that of the inverse bath temperature j3 in equilibrium. 
We now compute the integral in Eq. (jSJ) by steepest descent. The saddle point z*, if it exists, 
must satisfy the above constraint and the saddle point equation Q = G(z;t,t w )\ z=z *ut w \, with 

*. = rfti q{k ^:\ s - gU(t, t w ). (16) 

2 J A (2ir) d 1 - zq(k,t,t w ) 

Restricting the analysis to 2 < d < 4 and using Eq. (fl5l) . one finds that G(z = — flo) = Q c < 
is finite, while G approaches infinity as z tends to 0m- Therefore now the saddle point equation 
admits a solution only for Q > Q c , and the integration path is shown in the inset of Fig. [21 
Instead, for Q < Q c , exploiting the analyticity of h in the neighborhood of the branch point 
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FIG. 1: The large deviation function C(Q;t,t w ) is plotted against Q, in d = 3, for t w = 10 and two 
values of x for a quench to T = (continuous lines), or T ~ T c /3 (dashed line). Circles represent Q c . In 
the inset the same quantity is plotted in equilibrium at T < T c . 



z = —fto, an analysis similar to the classical one of [l2[ shows that the steepest descent route 
deforms into a cusp whose peak is sticked in z = —/3q (see inset of Fig. (j2J)). With this saddle 
point structure, finally one obtains 



2nP(Q; t, t w ) = exp[Vh(Q; z steep ; t, t w )\ 



^steep 



z*(Q;t,t w ) for Q > Q c 
-/3 for Q>Q C 



(17) 



The heat large- deviation function C(Q;t,t w ) = h(Q; z steep ;t,t w ) is plotted in Fig. ([T]). Due to 
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FIG. 2: The quantity q(k,t,t w ) is plotted in d = 3 against y (for large t w ) and for two different choices 
of x = t/t w . In the inset the analyticity region of h (Eq.©) and the steepest descent paths are plotted 
(dashed red). 



the sticking of z steep , it consists of two parts. For Q < Q c it is linear. For Q > Q c it grows to a 
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maximum and then falls off (again linearly for large Q since z* — > 0m)- The two branches merge 
at Q c with a discontinuity in d 3 C/dQ^ 3 . 

A singular behaviour qualitatively similar to that described above has been observed in the 
numerical simulations of the quenching dynamics of a disordered model for glassy systems j3] • This 
may suggest a certain generality of the phenomenon in aging system and a possible common origin. 
Non-analytical large deviation functions have been also found in a stochastic dissipative model for 



a single particle 17J], in simple non-equilibrium systems coupled with two reservoirs 18|, Il9[, and 



in diffusive mode 
behaviour in 



s in the continuum limit 



20 



2 If . always in stationary conditions. The singular 



17H19| has been related to the occurring of rare and very large fluctuations in the 
initial distribution. 

In the non-equilibrium setting considered here, the singular behavior of the distribution is 
related to the tying of the saddle point solution to the analyticity edge. This mechanism is 
mathematically similar to the one occurring in the equilibrium phase-transition. In that context, 
the zero wave vector fluctuations develop a macroscopic variance j3, through a mechanism 
reminiscent of the Bose-Einstein condensation. A similar phenomenon is dynamically produced 
here in the realm of fluctuating quantities: when a large amount VQ < Q c = VQ C of heat 
is released, a macroscopic fraction Q — Q c is provided by the k = mode. This is a novel 
condensation mechanism for non-equilibrium fluctuations. 

The large deviation function exhibits remarkable symmetry properties in the limit t w — > oo 
with Q fixed. It can be shown that in this limit the first term (i.e. — zQ) in Eq. ([9]) is dominant, 
implying that the above limit amounts to test the behavior of the tails of the heat probability 
distribution. In this regime one finds an expression where only Q = Q/qt yp and x = t/t w appear 

£(Q;t,t w ) = C(Q;x), (18) 

with C(Q;x) = qt yp /3oQ or C(Q; x) = qtypfiuQ for Q < or Q > 0, respectively. This shows that 
in this process the same scaling symmetry, which holds for average quantities, underlies also the 
behavior of fluctuations. As a consequence, C(Q;t,t w ) takes the simple form of Eq. ( 1181) when 
its arguments are measured in units of their reference value at the current age t w of the system. 
Moreover, using the expression of C(Q;x) one finds the asymmetry function 

£{Q; t, t w ) - C(-Q; t, t w ) = -{(3 M - /3 )Q. (19) 

Plugging this result into Eq. (HM one recovers a relation formally identical to Eq. ([T|). Notice 
however that the physical context is quite different: Eq. (CD) holds for t — t w large, while the 
validity of Eq. ( 1T4|) requires V large. Apart from this difference, Eq. ( Tl9l) shows that, by virtue 
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of a scaling symmetry, a fluctuation relation like (CQ) may be obeyed also in systems that are 
not at stationarity, but are slowly relaxing and aging. To the best of our knowledge, this is the 
first analytical result showing this in a classical model of statistical mechanics with a non-trivial 
equilibrium phase diagram. The quantities /3q,/3m represent the origin of the cuts of h, and in 
close analogy to the equilibrium case can be regarded as self-generated internal temperatures. 
According to Eq. ffl9|) these regulate large heat fluxes. Recalling that — (3q = = 0) and 

Pm = q~ Y (k = fcju), such temperatures can be naturally associated to the ordering modes releasing 

— * ~~ * — 1 /2 

energy at k = 0, and to those absorbing heat at a finite wave vector kyi = Vm^w (see Fig. [2]). 
Notice that, for t w — y oo and fixed x, /3 and f3 M decrease to zero as q^~y P = t]iT ■ Interestingly 
enough, this is the same behavior observed for the so called effective temperature e ff, defined in 



terms of the ratio between the response and the correlation functions, in the present model [15|. 
However, the relation between the quantities /3q, (3m entering P(Q;t,t w ) and e ff remains to be 
fully clarified. 

Finally, we briefly discuss the modifications introduced to the present picture by a quench to a 
finite temperature. It has been shown 15| that, in this case, the order parameter can be split into 
two statistically independent fields ip = cr + if), where cr and t/> are, respectively, an ordering and 
a thermal fluctuation component. In the scalar case (N = 1), these two terms are associated to 
the slow aging process caused by the displacement of interfaces and to the fast spin fluctuations 
with equilibrium character inside the bulk of the domains. This additive property amounts to the 
splitting C = + C^> of the correlation entering our calculations. In 15J it is shown that C^> 
is the equilibrium correlation at the quench temperature T, while behaves as in a quench to 
T = 0, apart from some trivial non-universal constants. Then, in place of Eq. fjSj) , one arrives at 

/ +200 + 2Q 
te V [HV + HW + W] t (2Q) 
-ioo+zo 

where h {a) and are given by Eq. flSJ) by setting the correlator C = C {a) or C = C w 
respectively, and hS a ^ is a function containing cross products C^CW. In the limit t w — y oo with 
fixed x, or alternatively with t — t w fixed, it is possible to show that the cross-term h^ a ^ can be 
neglected. Hence the heat probability results as the convolution 

/ioo+zo . 
dze v[hW+hW] (21) 
-ioo+ZQ 

of the fast and slow degree distributions. has the properties discussed insofar for the quench to 
T = while is the equilibrium distribution at the temperature T. Notice that, in the regime 
t w — y oo with t — t w fixed, is negligible and one remains with the equilibrium distribution 
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alone P = Pw. With these behaviors, it can be shown that the saddle point structure described 
above is not changed, except for a shift of the branch points at @o,/3m Then a singularity in the 
large deviation function occurs at a temperature dependent Q C {T) < 0. As it can be seen in Fig. 
[U the convolution with the equilibrium part, produces a broadening of £ particularly for large T 
and/or small x. This convolution structure, shown here for the first time, is expected to be very 
general in aging systems where a wide separation of time scales occurs, and also appropriate for 
other fluctuating quantities, beside Q. We notice that an analogue property is not expected in 
critical quenches at T = T c where the additivity ip = a + ifi is not obeyed: The composition of 
equilibrium and off-equilibrium fluctuations in this case remains an interesting issue to be clarified. 

By summarizing, we have computed the exact asymptotic probability distribution of the heat 
exchanged by a quenched ferromagnet described by the large- N model. A rich scaling structure 
emerges where heat, released by the small wave vector ordering modes, flows to components with 
finite wave vectors. The heat large deviation function shows a non-differentiable behavior with a 
singular threshold Q c signalling the onset of fluctuations condensation at zero wave vector. Heat 
currents are governed by a fluctuation relation analogue to the one obeyed in stationary systems 
in contact with two baths, but here with two self-generated temperatures fto, Pm- It is a challenge 
to establish to what extent the scenario above outlined is generic and holds also for systems with 
finite N. 

Acknowledgments 

We acknowledge fruitful discussions with Felix Ritort, Sergio Ciliberto and Leticia F. Cuglian- 
dolo. GG acknowledges support by PRIN 2009SKNEWA. 



[1] H. Touchette and R J. Harris, Large deviation approach to nonequilibrium systems, in R. Klages, W. 

Just, C. Jarzynski (eds), Nonequilibrium Statistical Physics of Small Systems: Fluctuation Relation 

and Beyond (Wiley- VCH, Weinheim, 2012). 
[2] L. Bertini et al, J. Stat. Phys. 107, 635 (2002); Phys. Rev. Lett. 87, 040601 (2001); J. Stat. Mech. 

P07014 (2007). T. Bodineau and B. Derrida, Phys. Rev. Lett. 92, 180601 (2004). J. Tailleur, J. 

Kurchan, and V. Lecomte, Phys. Rev. Lett. 99, 150602 (2007). 
[3] D.J. Evans et al., Phys. Rev. Lett. 71, 2401 (1993); G. Gallavotti, E.G.D. Cohen, J. Stat. Phys. 80, 

931 (1995). 



10 



[4] For a recent review see U. Seifert arXiv:1205.4176vl [cond-mat.stat-mech]. 
[5] C. Jarzynski and D.K. Wojcik, Phys. Rev. Lett. 92, 230602 (2004). 
[6] T. Bodineau and B. Derrida, C.R. Physique 8, 540 (2007). 

[7] A. Crisanti and F. Ritort, Europhys. Lett. 66, 253 (2004); F. Ritort, J. Phys. Chem. B 108, 6893 
(2004). 

[8] A. Imparato et al, Phys Rev. E 76, 050101(R) (2007). 

[9] J.R. Gomez-Solano et al, Phys. Rev. Lett. 106, 200602 (2011). 
[10] L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993); J. Phys. A 27, 5749 (1994). 

L.F. Cugliandolo, J. Phys. A 44, 483001 (2011). 
[11] A.J. Bray, Adv. Phys. 43, 357 (1994). 
[12] T.H. Berlin, and M. Kac, Phys. Rev. 86, 821 (1952). 

[13] C. Castellano, F. Corberi, and M. Zannetti, Phys. Rev. E 56, 4973 (1997); N.Fusco and M.Zannetti, 

Phys. Rev. E 66, 066113 (2002). 
[14] A. Coniglio and M. Zannetti, Europhys. Lett. 10, 575 (1989). 
[15] F. Corberi, E. Lippiello, and M. Zannetti, Phys. Rev. E 65, 046136 (2002). 
[16] R.J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic, London, 1982). 
[17] J. Farago, J. Stat. Phys. 107, 781 (2002); Physica A 331, 69 (2004). 
[18] P. Visco, J. Stat. Mech P06006 (2006). 

[19] R.J. Harris, A. Rakos, and CM. Schiitz, Europhys. Lett. 75 227 (2006). 
[20] G. Bunin, Y. Kafri, and D. Podolsky, arXiv:1207.6117 [cond-mat]. 
[21] L. Bertini et al, J. Stat. Mech. L11001, (2010). 



